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ABSTRACT 


A method is devised for the transient-state thermal 
stress analysis of two pipes, joined butt to butt and 
subjected to rapid or sudden change of internal fluid 
temperature. Although it 1s assumed that there is 
symmetry about the common pipe axis, the properties of 
the materials as well as the thickness of each pipe may 
be different. Then, given a specified time-temperature 
relationship for the internal fluid, over a specified 
problem time, the maximum stress encountered may be 
obtained. A digital computer program is appended for 


the solution of such problems. 
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inner radius of pipe 

outer radius of pipe 

specific heat 

plate flexural rigidity = Fh>/12 (1-1~} 
Young's modulus of elasticity 

sum of errors for least squares analysis 
pipe thickness 

coefficient of surface heat transfer 
thermal conduetivity 

bending moment (per unit circumference) 


dislocation bending moment (per unit 
circumference) at pipe junction 


total bending moment (per unit circumference) 
at pipe junction = ee 


bending moment (per unit circumference) 
remote from pipe junction 


number of data points in thermal transient 
analysis 


radial coordinate with origin at center 
of pipe 


a+b 
2 


time 
temperature 


radial displacement, remote from junction, 
positive outward 


volume, also radial force (per unit 
circumference) 


radial force (per unit circumference) at 
Pipe jume eon 


aXlal displacement, measured positive 
away from pipe junction 


radial distance, meaSured positive outward 
from mid-surface of pipe 


axial coordinate with origin at pipe junction 


Greek Letters 


coefficient of linear thermal expansion 
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J r“h? 

ene 
Fourier modulus = CoV 


unit elongation in radial direction 
unit elongation in circumferential direction 
unit elongation in axial direction 


radially outward deflection due to dis- 
location stress only 


tip slope = dn/dz 
Poisson's ratio 
weight density of pipe material 


thermal (radial, circumferential, axial) 
stress remote form pipe junction 


dislocation (circumferential, axial) stress 


total (radial, circumferential, axial) stress 
principal stresses 
maximum significant stress 


shearing stress 
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LT; INTRODUCTIGH 


In the design and analysis of certain piping of critical 
importance, such as nuclear power piping, approved methods, 
such as those in the USAS B31.7 Code, call for analysis of 
transient thermal stresses. The most important cause of 
such transient thermal stresses is variation of the 
temperature of the fluid flowing in the piping. It is 
conventional to regard the fluid as having infinite heat 
capacity and to presume that the fluid flows through the 
Piping with sufficient rapidity that all parts of the piping 
component or assembly in question are simultaneously 
exposed to the same interior ambient temperature. At any 
rate, an analysis based upon such simplifying assumptions 
may provide the engineering answer that is required or may 
form the basis for a more precise evaluation. 

The geometry of piping components and assemblies is 
generally rather complicated and one method of analysis is 
to study a greatly simplified geometry and then to apply 
"stress indices" which are tabulated constants intended to 
reflect the greater severity of the thermal stresses in the 
actual object as compared to that having the idealized 
geometry. One such simplified geometry of importance is 
that of two semi-infinite pipes, having a common axis, and 
joined together by welding or brazing at the origin of an 


axial coordinate system. If the two pipes have identical 


dimensions and properties, determination of the transient 
thermal stresses is a rather simple problem: there is no 
variation along the axial length. However, if the two 
pipes do not match in thickness or in material properties, 
so-called dislocation stresses will be set up near the 
junction and the task of determining the significant 
transient thermal stresses becomes considerably more 
difficult. It may be noted that the analysis to be given 
herein also contemplates a mismatch in pipe diameter, but 
this is not important since the theory employed will be 
limited to the case of thin-wall pipes. 

The problems of interest arise upon considering the 
ambient fluid temperature, to which the interior of the 
piping is convectively exposed, aS varying with respect 
to time in a specified fashion. In the development which 
follows, the general approach is capable of dealing with 
any temperature-time relation. The implementation presented 
is in the form of an immediately useful digital computer 
program for the analysis of these problems. In addition, 
there is specialization to temperature-time relations 
which postulate a long period at constant temperature 
followed by a sudden change in temperature or by a rapid 
linear variation in temperature which is in turn followed 
by maintenance of a new steady temperature. 

The analysis of the thermal stresses in any case con- 


sists of the following steps: 
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(1) determination of the temperature distribution 
in the piping 


(2) determination of the thermal stresses in the piping 
at a considerable distance from the junction 


(3) determination of the dislocation reactions, i1.e., 
the radial force and bending moment (per unit 
circumference) exerted by one of the semi-infinite 
pipes upon the other, and the calculation of 
stresses resulting therefrom 

(4) combining these stresses in an appropriate manner 
to find what could be called the "significant 
stress" and locating the point in the piping and 
the point in time for which this significant stress 
is greatest 

The task is greatly simplified because the problem at 

hand is one possessing symmetry about the common pipe axis. 
It is assumed that the fluid temperature is the same func- 
tion of time at all points adjacent to the inner wall of 
the pipe and that the coefficient of convective heat 
transfer is likewise the same at all such points. It is 
also assumed that the exterior surfaces of both pipes are 
perfectly insulated and that the following properties of 
both pipes are given and classified into two categories: 
first, thermal conductivity, specific heat, and density 
which may vary with temperature; and second, Young's 
modulus of elasticity, Poisson's ratio, and the coefficient 
of linear thermal expansion which must remain constant. 
Under certain conditons the properties in the second cate- 
gory might be permitted to vary, but this has not been 


investigated in this thesis. It is also necessary to know 


pipe wall thickness and outside diameter. 


ie 


The task of determining the transient temperature 
quistribution in each Cf Wene ewespipes iISesampieeeiued by 
the assumption that there is no axial flow of heat even 
though adjacent points in the two pipes may not be at the 
Same temperature. In other words, an error is introduced 
by neglecting axial flow of heat near the junction but the 
problem is too difficult to handle without this assumption. 
There is no evident way to estimate the effect of the 
assumption but the effect should be small for thin wall 
pipes to which this entire development is aimed. 

With this assumption, the determination is relatively 
straightforward. If the properties of the materials were 
not themselves dependent upon temperature, there might be 
some incentive to seek analytic solutions. However, for 
the sake of generality, a numerical method has been used. 
This consists of slight modifications to a digital computer 
program, oEriyr, developed by JOnn Ev earoeck, Frolressor OF 
Mechanical Engineering, Naval Postgraduate School, Monterey, 
California. This program is based upon the general algo- 
rithm described in Reference 5, but including several 
auxiliary features; see Appendix B. 

The output of BETTY] is a tabulation, at each instant 
in time, of the temperature at equally spaced radii through 
the pipe wall. In the present application this provides 
too much data to handle so that the next step is to pass a 
polynomial approximation of chosen degree (usually of third 


or fourth degree) through the data generated by BETTY1, and 
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retain the coefficients of this polynomial for each instant 
of time. This also provides for more Convenienkt procesemag 
at later stages of the analysis. 

Next, using this information, the stress distribution 
and radial displacement is calculated for each pipe at a 
position remote from the junction. This is a simple and 
standard sort of calculation which is simplified by virtue 
of the temperature data being given in polynomial form. 
Here, and elsewhere in what follows such calculations must 
be made for each pipe at each instant in time. 

The next step is the determination of the dislocation 
reactions, 1.e., the radial force and bending moment per 
unit circumference at the pipe junction. This calculation 
1s based upon the fact that continuity of radial deforma- 
tion and its first derivative with respect to the axial 
coordinate must be maintained at the junction. It leads to 
a pair of simultaneous equations in which the results of 
the calculation described in the preceding paragraph enter. 
Having the dislocation reactions, the corresponding disloca- 
tion stresses may be determined. These all contain a nega- 
tive exponential term and die out away from the junction. 
However, they are not necessarily maximum right at the 
Junction, so that some care is required later to be sure 
that the worst condition is not overlooked. The actual 
stress components at any point are obtained by superposing 
the thermal stresses remote from the junction and the dis- 


location stresses near the junction. 
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Finally it is necessary to determine the "significant 
stress" or measure of stress severity. In conformity with 
the stress criteria used by the most applicable codes, the 
maximum shearing stress criterion has been employed here. 
At any cross section, at any radius, and at any time, this 
maximum shearing stress! may be calculated without great 
difficulty. Because the variation with respect to radius 
occurs only in the non-dislocation stresses (except for a 
linear variation) and because the variation with respect to 
axial position occurs only in the dislocation stresses, the 
problem of determining a maximum shearing stress is 
simplified. 

In brief, the analysis and its implementation in the 
form of a digital computer program called PAM1 perform all 
the above steps and emerge with a statement of the magnitude 
of the "most significant stress," its location, and the 
titie atowhich itjsoceurs. 

The effort in developing this thesis has been in under- 
standing and relating the several steps as described above, 
in writing sub-programs to perform the operations in each 
of the steps as outlined above, making use of programs 
(such as BETTY] and a least squares curve fitting routine) 
that appeared to be of immediate use, debugging and proving 


each portion, knitting the separate parts together into a 


ee is conventional to calculate and employ a quantity 
known as the "significant stress" or the "Stress intensity" 
which is equal to twice the maximum shearing stress. This 
practice is adhered to in what follows herein. 
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complete whole, and finally debugging and proving the final 
program. 

The steps of the theoretical analysis and the methods 
devised to employ them are described in Sections II and III 
hereof. The digital computer program is described, in 
general terms, in Section IV. Instructions for the use 
of the program are given in Appendix A. The program list- 
ing is given in Appendix B. 

There was not sufficient time after developing the 
program to make a systematic use of it to investigate the 
effect of certain material or geometrical parameters upon 
the value of the significant stress. It is recommended 


that this become the point of departure for a future thesis. 


YS 


[ot eee Olen, 


When the wall of a pipe is nonuniformly heated as for 
instance when it is subjected to a rapid or sudden change 
of internal fluid temperature, its elements do not expand 
uniformly and mutual interference sets up thermal stresses. 
It is assumed that the distribution of the temperature does 
not vary along the axis of the pipe and that there is 
complete symmetry about the axis. A ring is cut from the 
pipe by two plane cross sections perpendicular to the axis 
and a unit distance apart. During thermal deformation such 
cross sections can be assumed to remain plane if taken 
sufficiently distant from the ends of the pipe. Therefore, 
the unit elongation in the direction of the axis is constant 
[7]. Then, the unit elongations in the perpendicular 


qimect lemoaware 


ee = dw. constant 
Fy az 
du 
= — il 
a ar (1) 
3 —_ UU 
£ Ve 


These elongations can be expressed as functions of the 
axial, radial, and circumferential stresses remote from the 
end of the pipe Orr ee and oe respectively), the thermal 


expansion, %, and the temperature, T, where T varies with 


Piemradiial di Stance, ©£,. only. 
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From Reference 6, p. 66, the elongations are 


O 
ee ee ee 
a S F ‘On, ‘ O,) ue 
oe r 
Ex = [7 |; (oO, ae O,) oor: (2) 
O 
ie LU 
= ee + + 
a = = (o,, o,.) aT 


hee A) dew teacammerease In Unit volume. From tunis 


(G oe } eee a (3) 


Zz, ie te E Z ie oe 





Substituting equations (2) into equation (3) it can be 


Shown that 














E u aTE 
= + == 
oF 1+u cy it ae. l- 2 
E u aTE 
= + — 
CO, Ta (e. ean T- 2 (4) 
ls u os 
= Sa es ee 
2 amare ees oi. Lao 


From Peference 7, p. 206, the equation for radial equi- 


imeem Of an element is 


Substituting equations (4) and (1) into equation (5) 


ea 1 du u 1l+y ' ae 
ar? r “dr 2 ip dr 
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which may be written in the form 


ee Ab 


a 
[ ae ae = ene ec (7) 





A [e 


Ga 
Glia 
Integration with respect to r gives 


d 1 +4 
ae =— + 
an (ru) ia Oe 2C 4x (8) 





A second integration gives 





i 
1 sl+hn il 
=< + + oe 
u =e d aTr dr + Cyr + C, ‘ (99) 
where Cy and C. are constants of integration. Cy and C. are 


found from the boundary conditions at the inner and outer 


surfaces of the pipe, namely 


a. a > ee (10) 


From equations (10) the constants of integration are eval- 





uated as 
1 ‘i a2 b 
C = ee Jee e 
- ae a = a’ : . 
(1 + w) (1 = 2) i 2 
Cy = i ne) fe aTr dr - be, (11) 
u b”" - a“ a 
An expression for Die 1s obtained by substituting x a ay Gia 
and Ey = u/r into the second equations (4) and then taking 


u from equation (9) and the values of the constants from 


equations (11). This gives 
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F 1 iG oo _ ze b + 
= = oo i Gprereme Clie ar, > a io Tr sais 
Tamer rob” -=sangEee 


— . 


el 





An expression for o, is obtained from the equation of 


t 


equilibrium (5), thus 





&, wcities 
- io Cir eis 
From this, 
vs Z 2 b j 
Oo, = 2 i ey Corer aa Sf aTtr-dv— areas 
ies l- 2 Zane 
ie a r (b - a’')a 


An expression for oF is obtained by substituting ce. = du/dr 
and Ee, = u/r into the first of equations (4) and then taking 
u from equation (9) and the values of the constants from 
equation (ll). This expression will contain the constant 
elongation, oe in the direction of the axis of the pipe. 

If it 1s assumed that the pipe can expand freely in the 


axial direction, the magnitude of Ee, can be calculated 


from the condition that the sum of the axial forces over 


any normal cross section 1s equal to zero. From this 
D b 
is — ore ar (14) 
b* - a a 
and 
b 
oe = = J Vatr rds cs (15) 
(1 Sap los a a 


Ihe, 


Substituting the values of the constants from equation (11), 


the radial displacement remote from the pipe junction is 





2 
a a 
U =—giren tol ) — f aTrdr + —~—~— f oaTrdr! (16) 
Z 1l- wu r eg 2 7 
a b -a a 
Thus the thermal stresses, oe Oy and CoG and Gadaal 


displacement, u, at a distance from the pipe junction may 
be determined by use of equations (12), (13), (15) and (16). 
It should be noted here that at a distance remote from the 
pipe junction there is no variation of the stresses or 
radial displacement with respect to z, the axial coordinate 
of the pipe. There exists here a bending moment per. unit 


circumference, denoted as M,, where 


b 


Me Cement ee clr (17) 
a 





Zz 2 


See Figure 1, case ae 

On the other hand, if the end of the pipe is supported 
or clamped, as for example when it is joined to another 
pipe, free expansion of the pipe is prevented and a bending 
moment per unit circumference, Mo and a radial force per 
unit circumference, View are set up at the pipe junction 
(Figure 1, case B). Because of the continuity with the 
joining pipe, the moment, stresses, and radial deflection 
at the pipe junction will be different from those at a 
distance from the end of the Pipe. @in other words, the 


adjoining pipe exerts a bending moment and radial force. 
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It should be noted here that the total bending moment at 
the pipe junction is M. - M., which will peewee wecemeed 
here as Mo: Thus, in solving for the dislocation reactions, 
although Mo is obtained directly, it is worth noting that 
that bending moment depends on both the moment at the pipe 
Junction as well as that remote from the junction. The 
total radial force is just that at the pipe junction, Vive 
Since there is no force remote from the junction (Figure 1). 
This bending moment and radial force cause a aecileeren 
and tip slope at the end of the pipe (z = 0), which are 
denoted by n and 6, respectively. The same is true, however 
for the other. Therefore, let ny and Ol, denote the deflec- 
tion and tip slope, respectively, at the end of the adjoin- 
ing pipe. Now if these pipes are to join up continuously, 
the total deflection for each pipe at the pipe junction 
must be the same. The total deflection consists of the 
deflection due to M and Ve as well as the radial displace- 
ment, u, remote from the junction. Equating these total 


deflections, thus 


ntue= ny + Uy (18) 


Also, if these pipes are to join continuously, the tip 
slopes at the ends of each pipe must be negative of each 


other. That is, 


A = = 6 (18a) 


Ze 





IY wh CASE A (Independent of 2Z) 





> CPE FO RAE A Oy Re a 
oe - —- wes oe rh Me I Eee hee EE AE, CEE - OPES warren ane ae NE PT et TS, OE ee ee ee 
( vm es —— ee ean — a 
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Figure 1. Combination of Cases A and B 
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From Reference 7, it is seen that 





-Bz 
n= = 5 [BM (-Sin 8z2+cos 82) + V_ cos 862} (19) 
O O 
Zo) 
where 
Fh? 
as 
12(1 - uw’) 
and 


‘ite 
a sche 


——)“h 

( 5 ) 
and where the algebraic signs have been changed to conform 
to the convention herein that n is positive outward. If 


the pipe junction is taken as the origin of the axial 


coordinate (i.e., z = 0), then the deflection at the junc- 
tion is 
|) = — [BM + Vo] (20) 
Zp PD 


Also from Reference 7, the tip slope is the first derivative 


of the deflection with respect to Zz, or 


dn _ aa 


dz 28°D 





[28M. cos £2 + Ve (cos 82 + Sin 82)] (21) 


Therefore the tip slope at the pipe junction (z = 0) is 


1 


28°D 





© 
I 
l 


[2BM. ess (225) 
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It should be noted that the radial displacement, u, 
given in equation (16) is a function of the radius, r. 

The value of u used in equation (18) was the average value 
of the radial displacement found at the inner and outer 
radii of the pipe. This seems to be the most reasonable 
and convenient of several different choices that could be 
made. 

Noting that Bg and D are functions of the material and 
the geometry of each pipe, their values may be obtained and 
substituted into equations (20) and (22). Also the values 
@£F uU and u,, May be obtained from equation (16). For 
mechanical equilibrium, we must have Mo= M'- M =", =8)) sen 


O oo 


and Ms = -V These relations, together with equations 


Oils, 
(18) and (18a) permit solving for the dislocation reactions 


M', M' 


O ei and i Hence, the dislocation stresses near the 


pipe junction can be obtained. 

From Reference 7, the axial stress due to the bending 
moment is 

es eee YM) (23) 


ZZ h> 


where y is the radial distance measured from midsurface, so 


felled a - WS x. Also from Reference 7 
Yo -BZ ' ' -BZ ° 
M(z) = [fe Sin Bz] + Mo [e (cos 8z+sin8Bz) ] (24) 
The maximum axial stress is obtained at the axial 
position where M(z) is a maximum. Since M(z) contains a 
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negative exponential term, it will eventually decrease to 
zero aS Z increases. But this does not necessarily mean 
that M(z) will be maximum at the origin. The method of 
Finding the point of maximum bending moment is shown in 
detail in the next section. 

From Reference 7 the circumferential stress due to the 


eas locatrvenereaceions 1S 


_ En(2) 
Ope + HOLS (25) 
m 
Where n(z) comes from equation (19) and ee =. ies li we) 


be noted here that there 1s no radial stress due to the 
dislocation reactions because of the assumption that the 
Pipe wall is thin. 

Summarizing now, there is a radial, circumferential, 
and axial stress remote from the junction as well as 
additional circumferential and axial stresses near the 
junction. Each of these is some function of the radius, 
so that for a given radius and given axial location these 
stresses may be superposed on each other and the total 
stresses found. Using equations (12), (13), (15), (23), 


~ 


mare 25>) these total StresSesmare detinea as 


ae a 
Yr 
Onn = 01. + O., (26) 
tc 
= + 
On oD O32 


iaeceieadwng Subseript fF indicates “total.” 


a 


The radial force V, also gives rise to a shearing 


stress, Presuming that this shearing stress is 


a ° 
eee 


distributed parabolically, it can be represented as 


_ 3V _ ,2y,¢ 
tmres — aN [1 (<*) ] (2970) 


where, from Reference 7 
VS a ie (cos Bz - sin Bz) - 28M. Sin 82] (28) 


Thus a general element is loaded as shown in Figure 2. 
(It should be noted that this element is also subjected to 


the circumferential stress, d,,). 





Figure 2 
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The three principal stresses are then 


Soi Oo oe 
O O io - oO D 
oS ee ie Tz Tyr 2 
So - i a eae (29) 
" 
Foe one, Onmy 7 (2 ee, 5 dae Z 
3 | 2 Trz 


The maximum significant stress at a point is defined as 


sig = Max {lo, - o,/,|0, - 93|,|9, - 9, ]}. (29a) 


It is necessary to find the point in space (i.e., values of 
meand 27) and the point in time for which cee achieves its 
maximum value, and to determine this value. This must be 
done for each of the two pipes individually since allowable 
values of O..4 may be different for the two pipes. How this 


1g 
Mmemaone 1S considered in detail aiaSectiemel il. 


a 


ee DE TOD 


Se iindecateds.1n SSCeHone, Ene preqranetEtauhas been 
incorporated into the program PAM] to obtain the transient 
temperature distribution. This is a program for dealing 
with the one dimensional thermal transient analysis of 
homogeneous slabs or tubes. Its use herein is for the 
thermal transient analysis of a homogeneous pipe in which 
the outer surface is perfectly insulated and the inner 
surface is convectively exposed to a fluid having a 
specified time-temperature relationship. Physical prop- 
erties of the pipe material (i.e. thermal conductivity, 
weight density, and specific heat) may be temperature 
dependent; the program is set up to consider property- 
temperature relationships approximated by a third degree 
polynomial. A similar variation 1S permitted in the 
surface heat transfer coefficient. 

The following is a brief summary of the general 


algorithm of the program [5]. As shown in Figure 3 a unit 


length of pipe, contained in a small dihedral angle, is 
considered to be divided into subvolumes, where the inner 
and outer subvolume have half the thickness of the others. 
The thermal properties of each subvolume are considered to 
be concentrated at its central nodal point, and heat is 
imagined to be conducted between nodal points through a 
network of fictitious heat-conducting rods of appropriate 
thermal conductance. In the transient state, however, 
heat is not only conducted to and from each nodal point in 
the lumped network, but additionally each nodal point 
experiences a change in internal energy in an amount which 
depends upon its temperature change during a given time 
increment, its specific heat, the total volume of material 
which it represents, and the density of the material. 
Assuming that the specific heat, c, and the thermal 
conductivity, K, are uniform, the heat conducted by each 
rod is -KAT, where AT is the temperature difference 
between adjacent nodal points. The total heat conducted 
in the finite time increment At is -KATAt. Now the change 
in internal energy of a given nodal point is covAT , where 
o is the material density, V is the subregion volume, and 
T is the temperature of the nodal point at the end of 
the finite time interval At. The heat balance can be 


written as 
=~LKATAC = "evar (2.0) 


Using an example of three nodal points, as shown here 


Zo 





C,0,V,K 


| 
| 
| 


Figure 4 


a specific example of equation (30) is 


K(T, = Ty) ate K(T. = T,) a 


SolVinigmeror ar the future temperature of the nodal point, 


Ale 
TT acacia eee = et la (31) 
1 2 3 Y i 
where 
2 KAt 
Y Cow 
and is known as the Fourier modulus. From this it can be 


seen that the future temperature of 1 depends on its 
present temperature and the present temperature of its 
adjacent nodal points. Equation (31) is the basis of the 
transient-state numerical method. It should be noted 

that the primary purpose of BETTY] was to find the quan- 
tities delta-tee-one and delta-tee-two, defined and required 


by the Nuclear Power Piping code, but in the present 


30 


application only the transient-state temperature distri- 
bution is of interest. 

The immediately applicable accessible output from 
BETTY] is a set of nodal point temperatures at each epoch 
of time. For convenience in what follows, it was decided 
to represent this temperature distribution by a polynomial 
in r, using least squares procedures. This results in an 
economy of computer storage and convenience in performing 
the integrations required in equations (9) and (11) 
(et.seq.), with minor loss of accuracy. A brief discussion 
of the curve fitting procedure follows. 

The deviation at a data point is defined as the differ- 
ence between the tabulated T-value and the T-value computed 
for a functional relationship. The least-squares criterion 


asks for a minimum of 
: Z 
ie (x. ) = ial (32) 


where N equals the number of data points, T' denotes a 
value computed at r= rs from a functional relationship, 
and T; denotes the corresponding tabulated value. This 
can be differentiated to find its minimum, which leads to 
equations that in many cases of practical interest are 
linear and, in principle at least, easy to solve. 

The form of the relationship must be chosen in advance. 
Bimce it 1s important to-find a relationship which wili@be 


easy tO Obtain and integrate, a polynomial of the form 


ot 


i 
a oa 33 
j (33) 


was selected. Here k is the degree of the polynomial and 
the Goeffiierenes Agr Ayr Aan vee A must be determined to 
Minimize G. 


Substituting equation (33) into equation (32) 
(34) 


For a minimum, the partial derivative of G with respect 


to each coefficient, an! (ne —0; 1, 2 eek k) must vanish. 
Thus 
N K 
0 = a = 2 [2a > a r? - ey ae | 
a — a a 
1=1 J=0 
K x 1+m . m 
=2[2 a.(2 rs ) el (35) 
43=0 J i=l i=l 
m= 0, l, 2 k 


This set of equations may be arranged as 


N Dae > ees a ~T 
a I ali O st 
Le. 2G cro... Erstt a Pe. 
ile 1 ili jae 
a ; (36) 
ie ryktl pyrite 2k ay LT rs 
ab L al 


a2 


where, in this last expression, £ means 1 EEE as 
straightforward to calculate the ie te and to solve 
the system to obtain the desired coefficients a, 

For the sudden or rapid temperature changes which 
appear to be of greatest interest, it seems likely that 
the degree of the polynomial need not exceed k = 4, al- 
though no verification of this conjecture has been under- 
taken. However, it should be noted here that even for 
some unusual time-temperature relationship such as a 
temperature which fluctuates with time, a polynomial of 
larger degree could still be used to represent this 
temperature distribution. 

Now, for each pipe and for each time epoch, the co- 
efficients of the polynomial representation of the 


temperature as a function of the radius are available. 


These polynomials can be integrated and the stresses 


remote from the pipe junction can be determined. Equa- 
peat (1.2 ) 
is il Ma Ea : 
— - (a—- aor Ear + je ecr deals|) fle? ) 
7 an r* a r* (b* - a“) a 


will serve as an example here. Given values of E, u, a, a, 
and b, and the particular radius, r, as well as the poly- 
nomial for the temperature distribution, the radial stress 
can be found. On a computer, the radius can be incremented 


from the 1nneprepeatus, a> tOutkhevouter radius, byotheres, 


Bo 


finding the radial stress at increments throughout the 
thickness of the pipe. A similar procedure is followed 
using equations (13) and (15), thereby finding the circum- 
ferential and axial stresses throughout the thickness of 
the pipe. 

It was noted earlier that the radial displacement, 
remote from the pipe junction, used in determining the 
total deflection would be an average value of the radial 
displacements found at the inner and outer surfaces of 


the pipe, i.e., 
u = = fu(a) + u(b) ] (37 


From this it can easily be shown that 


b 
i (ath) f-—we -L ZOE UE f aTrdr] (38) 
2 Z b2 Z 
nc 
where E. is found. from equatwon (L4)eeinus cias “radram 


displacement can be determined for each pipe by evaluating 
the integral of the polynomial for the temperature dis- 
tribution at the inner and outer surfaces of the pipe. 

The next step is to find the bending moment, My and 
the radial force V, at the junction of the two pipes. This 
procedure was outlined briefly in the preceding section, 
but will be shown in some detail here. Substituting equa- 


tion (20) and the values obtained for u and u from 


Ales 


equaeron (38), into equation (18) "thus 
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1 ' il ' 
[6M + V_. i) u = [B.M eV lee 
28°D O O 28°D, Leul O il 


where as before, subscript 1 is used to distinguish between 


the two pipes. Now substituting Mo = M toe eect rT! 











Mol = M, e Moy! and rearranging 
M M 
il oo 00 
eo — ier —sealllemaaial 5 —ele 
Zoe D 2B,D, Zp 0) 2B,D, ZipaeD 2B,D, 


Substituting equations (22) into equation (18a) 


1 ' il i 
[28M + V_] = —= [28.M + Vj 
28°D O O 2B°D, LoL O 





then substituting as above and rearranging 








M M 
1 1 7 it 2 20] 
[—- + ——] M +[ + ) Vo. = -[—— + ~-—~=] (40) 
BD ByD, Oo '3g4p 285D, O BD B,D, 


Equations (39) and (40) are easily solved for the mechanical 
reactions My and Me, after which the appropriate dislocation 


heagectlion moments 


may be determined. 


Now that the dislocation reactions, Mo Moye 


have been determined, the stresses which they produce in 


and! Va 
O 


SS 


each pipe can be computed. These stresses may then be 
superposed on the thermal stresses already obtained and 
shown in equations (12), (13) and (15), obtaining the total 
stresses as defined in equation (26). Next computing the 
shearing stress, equation (27), and the principal stresses, 
equation (29), the significant stress may be determined. 

It is worth repeating here that this significant stress 

is the largest difference between any pair of principal 
stresses. Hence, it is a function of the radial coordinate, 
as well as the axial coordinate of the pipe. Also it may 
be useful to know which pair of principal stresses produce 
this maximum. 

The computer program takes the absolute value of the 
difference of each pair of the three principal stresses, 
and then selects the largest of these as the maximum 
Significant stress. This is done at increments of the 
radius throughout the thickness of each pipe. It has been 
noted earlier that the radial stress is zero at the inner 
and outer surfaces of the pipe. Also the circumferential 
and axial stresses remote from the pipe junction are maximum 
at either the inner or outer surface for a large class of 
fluid temperature variations. Therefore, it seems reason- 
able that in taking differences, the largest difference 
might be at either inner or outer surface. In fact, in 
certain important examples of specified time-temperature 
relationships (i.e. sudden change in internal temperature) , 


it has been observed empirically by the writer that this 
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maximum does occur at either the inner or outer surface, 
but admittedly this may not always be the case. [In any 
case, finding the radial position at which the significant 
stress is maximum is simplified by merely finding the sig- 
nificant stress at given increments of radius throughout 
the pipe thickness. 

The axial position at which the significant stress is 
maximum is that position at which the dislocation stresses 
are maximum. As noted earlier, these stresses are maximum 
at the point of maximum bending moment which may occur at 
the junction (z = 0) or at a mathematical extremum of the 
bending moment. Equation (24) gives the moment as a func- 
Gren, Of 2, £8, Mo and Vo: But now 8, Mo and Ve are known. 
Therefore, taking the derivative of the moment with respect 
to z and setting it equal to zero, the points of extremal 
moment can be determined. 


V 
raz) = a fom seen + M en awe Bz + sin®Bz) J (24) 


Setting the derivative equal to zero, 


O= NES maior Bz ~- sin £z)]J- em te P? (2 Sin Bz) ] 


Rearranging to a form readily usable on the computer, and 
setving for Z 


> =a a (41) 
B ae 28M, + V_ 


om 


This method is employed in the computer program. In 
this manner, the first (and only physically meaningful) 
point of extremal moment can be found for each pipe, being 
Sure to use the appropriate value of g. It must be 
remembered here, that the value of the dislocation moment 
atthe pipe junetionm Moe has already been determined. In 
examining certain important temperature-time histories, 
such as those corresponding to sudden or rapid changes in 
internal fluid temperature, it has been observed empirically 
that the maximum bending moment for each pipe occurs at 
the pipe junction. From these examples, it was noted 
that the magnitude of the moment, found by substituting 
the value of z obtained from equation (41) into equation 
(24), was always less than the magnitude of the moment at 
the pipe junction. In fact, in all the examples studied 
in this thesis investigation, not only is this moment 
at the junction a maximum, but the moment decreases 
rapidly away from the junction, oscillating above and below 
a value of zero. The program uses either the end moment 
or the moment at the axial position given by equation (41), 
whichever is larger, in computing the dislocation stresses. 

If the temperature distribution were such that the 
maximum significant stress did not occur at the inner or 
outer surfaces, the user of program PAM1 still has the 
option of wholesale determinations at different radii from 
which a maximum may be found. If the maximum significant 


stress did not occur at the end (z = 0), the preceding 
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determination, which is based upon maximizing M without 
considering n or V would probably not give the axial 
position of greatest significant stress. Because the cases 
which seemed to be of greatest interest lead to maximum 
Significant stress occuring at either the inner or outer 
surface, the program PAM] in its present form does not 
include the capability of making Jepasn ea at the 
axial positions other than for (z = 0) and the value of 

Z given in equation (41). Accordingly, for cases in which 
the maximum significant stress occurs between the inner 

and outer radii, the program does not presently provide a 
capability of finding the critical spot and determining 

the true maximum significant stress. However, there should 
be no great difficulty in modifying the program to provide 
this capability if desired. 

Once the maximum significant stress is found for a 
given time, the same procedure can be followed for the next 
time increment. It may be worth repeating that two values 
of maximum significant stress are determined at each given 
time epoch-one for each pipe. The temperature distribution 
is known for that time increment, therefore, a polynomial 
can be found to represent the temperature, and then that 
polynomial can be used to find the next significant stress. 
This method is employed in the computer program making it 
possible to find a maximum significant stress over any time 


period. 


39 


Lieb LG lTAtime COMPUTER PROGRAM 


This section will serve as a summary for all that has 
preceded it in the sense that now all the theoretical 
techniques and the methods devised to carry out these 
techniques have been incorporated into the program PAM1 
which is immediately usable for the solution of the partic- 
ular problem at hand. This program is written in Fortran IV 
language and employs only normal Fortran functions. 

Appendices A and B will give a complete listing of the 
program, as well as instructions for its use. The writer's 
purpose in this section is to outline the program (a flow 
chart in words perhaps) thereby describing the various 
Sub-programs employed and woven together to form the final 
program 

As indicated in the preceding section, the whole idea 
of applying the theoretical analysis devised herein to a 
transient-state problem rests on the premise that a tran- 
Sient temperature distribution could be obtained, and it was 
shown that by using a numerical method, such a distribution 
can be obtained. In view of the fact that this temperature 
had to be integrated in order to find the stresses which 
it produces, the temperatures and the corresponding radii 
found from BETTY1, were used as data points; then a poly- 


nomial was fitted to these points. The program employs the 
least squares method and Gauss elimination to find the 


coefficients of the polynomials. 
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PAM1 reads data for one pipe and solves for the co- 
efficients of the polynomial for every instant in time at 
which a temperature distribution has been determined. 

This continues until the problem-end-time- has been reached, 
with the program storing each of the coefficients thus far 
obtained. Next, the program reads data for the joining 
pipe and solves for the coefficients of a similar polynomial. 
Now, however, at each time interval the program recalls 

the coefficients for the other pipe at the corresponding 
time. The program continues to do this until the compu- 
tation terminates. For each time increment there is a 
polynomial for the temperature distribution for each pipe; 
hence, the stresses can be determined for every epoch of 
time. 

As shown in Section II, these stresses are determined 
by integrating the polynomial, which is a function of 
radius, multiplied by the radius (i.e. fTrdr). This process 
1s carried out in Subroutine PINT. Subroutine PVAL is then 
used to evaluate the integral at any given radius, as well 
as at the inner and outer surfaces of the pipe. [In cal- 
culating the circumferential and axial stresses, the 
temperature at a given radius is required. Subrouting TVAL 
evaluates only the polynomial for the temperature distribu- 


tion in order to obtain the temperature at a given radius. 


re ates ae ; 

This 1S an upper limit, introduced as an input to 
PAM1L. The computation terminates when time has been incre- 
mented to the value of problem-end-time. 
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The same techniques are employed in finding the radial 
displacement remote from the pipe junction. Subroutines 
PINTL, PVAL1, and TVAL1 perform the same procedures in 
finding the stresses and radial displacement for the 
adjoHning pipe’ 

The program then solves equations (39) and (40) for 
the unknowns Mo and Vor the total bending moment and radial 


force at the pipe junction, and the corresponding disloca- 


J 
Ole 


At this point, equation (41) is used to determine the point 


tion reactions (Mis ve) and (M Vv.) are next determined. 
of maximum bending moment. The moment at this point is 
then computed and compared to the moment at the pipe 
Junction, and the one of greater magnitude is used in 
determining the dislocation stresses. The total stresses 
are then determined and the maximum significant stress 

as defined in Section II is obtained. 

Each of the stresses obtained thus far depends upon 
the radius, so that the program can find these stresses 
for each pipe at any radius throughout the thickness of 
the pipe. Since it has been this writer's experience 
that the maximum significant stress occurs at either the 
inner or outer surface, the program offers the option of 
either finding only those stresses at the inner or outer 
radii or all stresses at given increments of radius. This 
depends on the wish of the user and is elected through 


the use of a data card. 
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In any case, once these stresses are determined for a 
given time the program continues to the next instant in 
time, calculates the coefficients of the polynomial for the 
second pipe, recalls the coefficients already determined 
for the first pipe, and performs all the evaluations 
described above. 

These procedures continue until the problem-end-time 
1s reached. After determining the significant stress for 
each epoch of time, the program calculates the largest of 
these, thereby determining the maximum significant stress 
encountered over the problem time for each pipe. PAM1 also 
furnishes the time and location of this maximum stress, 
as well as indicating the nature of the two principal 
stresses responsible for it. It should be repeated that 
the output includes two such combinations of results - one 


for each pipe. 
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APPENDIX A 


Instructions for use of Program PAM1 


This program is designed for the IBM series 360 com- 
puter; however, there is very little machine dependence. 
Testing and execution thus far has been on an 360/67 con- 
figuration. All quantities beginning with the alphabetic 
letters I, J, L, M, and N are integer variables, used for 
distinguishing specific locations in the pipe wall, used 
as counters for storing purposes, and used as limits for 
DO-loops. All other quantities are continuous (REAL *8, 
i.e., double precision) variables. Many of these variables 
differ in appearance only by the number (1) placed at the 
end of the variable. This is to distinguish related 
variables for each of the joining pipes. For instance, 
ELAS is used to denote the modulus of elasticity for one 
pipe, while ELAS] denotes the modulus of elasticity for 
the joining pipe. 

The program determines temperatures at the centers of 
each of the N layers into which the pipe wall is divided. 
The arrangement of these layers is one half-layer, (N-1) 
full layers, and a final half-layer. An example for N=5 


layers is illustrated in Figure A-l. 


44 





Figure A-1 


Note that the centers are numbered from C. to C There- 


Jl 


fore, in the program the given value of N is increased by 


aes 


unity so as to permit numbering the centers as illustrated 
in the figure above. One important consequence of this is 
that in fitting the polynomial to the temperature distri- 
bution, there are N+l temperatures and N+l radii used as 
data points. The number of layers should be between four 
and twenty, inclusive, but if greater than ten, it should 
be an even number. Otherwise the program automatically 
changes N to conform to these conditions. 

The temperatures, coefficients of the polynomials, and 
related stresses are calculated at discrete epochs, uni- 
formly separated by time-interval DT. The given value of 
DT, as selected by the user is reduced if necessary to 
assure convergence of the computation. As this calculated 
value may be different for each of the connected pipes, 
care must be taken to see that the smaller value is used. 
Because of the sequence of operations involved, it may 
happen that temperature evaluations have been completed 
for one pipe before analysis of the parameters relating 


to the second pipe has revealed that a shorter time 
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interval should have been used. In this case an appropriate 
message is printed and the computation ceases; it is then 
necessary for the user to make input modifications as 
advised and to restart the program. However, if this 
case does not arise, the computation procedes to completion 
as desired. Another method of insSuring that the same value 
of DT is used for both pipes is to remove temporarily the 
card DT = BB which precedes statement 31. This enforces 
the input value of DT. However, if this is done, one 
should be alert to the possibility of instability or non- 
convergence in the calculation. 

The data deck is arranged in two parts. The first 
part contains material and geometrical properties which 
are constant for each pipe. This part consists of two 
data cards. The second part of the data deck contains the 
numbers of layers into which the ee is divided, time 
interval, problem-end-time, data for those properties of 
the material and fluid which may vary with temperature, 
as well as information regarding the internal fluid time- 
temperature relationships and initial pipe temperature. 
The degree of the polynomial desired should also be 
selected at this time. Although there has not been 
sufficient experimental use of the program to warrant 
making definite statements, tentatively it may be recom- 
mended that the degree of the polynomial should not 
exceed 0.6 N (where N is the number of layers). Unless 


the fluid temperature fluctuates with time there seems to 
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be no reason to employ a polynomial of higher than fourth 
degree. Another important input is the option parameter 
NSTAR which determines if the stresses will be calculated 
at only the inner and outer surfaces or at given increments 
of radius throughout the wall of the pipe. Presently 
evaluations are made for r = a + n(b-a)/100; n=0, 1,...,100. 
However the interval may be modified by changing cards 900 
and 901. 

Detailed instructions for the preparation of each data 


Gara follow. 


meant I Card l 
Column Format DESeGEPOet Cimon Vomue 
ae) Bk. 0 Young's modulus of elasticity, 
psi 
20 fer” coefficient of linear 
: saul 
thermal expansion, °F 
21-30 ror 0 Poisson's ratio 
31-40 Eo 0 inner radius of pipe, inches 
41-50 Ee ono outer radius of pipe, inches 
Cancer 


This card is punched exactly as Card 1 using parameters 


now for the connecting pipe. 


In the following, note that all integer inputs must be 


right-justified in their respective fields. 
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Part II 


Column 


Zee 30 
31-40 


41-50 


Format 


eee 


3 


i 


tak 


az 


T2 


aS 


Ee O74 


F10.4 
Oe 


F10.4 


Card 


Cara 
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3 


Description or Value 


degree of polynomial desired 


1 or 2. Number one (1) is 
used to denote first pipe; 
number two (2) is used to 
denote second pipe. 


1 or 2. Option parameter one 
(1) used to calculate 
stresses at inner and outer 
surface only, while option 
parameter two (2) calculates 
stresses at increments of 
radius throughout pipe wall. 


Ll. {A digit in this positirer- 
cards 5 "through lil merely 
sets up appropriate internal 
branching) 


number of layers into which 
pipe thickness is divided 


1. Number one (1) must be 
punched here. Program BETTY1 
had capabilities to print 
output only on certain time 


intervals. Program PAM1 dces 
lal@ie: 


blank 


outside diameter of pipe, 
inches 


pipe thickness, inches 
time interval, seconds 


maximum problem time, seconds 


Co.) uma 


i! 


i= 2 0 


Pel 30 
31-40 
41-50 


E20 


Zr 30 
31-40 
41-50 


Wl 2 () 


SoU 
S40 
41-50 


Format 


i 


Poe 4 


F10.4 


dil 
F10.4 
F10.4 


Fiv0.-4 
F10.4 


Tl 
F10.4 
EO 


ELO: 
| all OF 


Dm 


Card 


Gard 


Card 
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5 


Description or Value 
2 


Surface heat transfer co- 
efficientyeeeu/hr = SO tue 
o 

EF 


Higher coefficients which 
give variation with respect 
to temperature according to 
the formula; Property = 

Bee oe EG: XX 7 ee eons) 
where X denotes temperature 
and F5 through F8 denote the 
four entries in columns ll 
through 50 respectively. If 
the higher coefficients are 
left blank, a value of zero 
is used. 


3 


Ene rMalk CONGUCETVity, eoeuy mr 
- ft - °F 


higher coefficients used as 
aepeeG asco 


4 
specific heat, Btu7lb = °F 


higher coefficients used as 
igee aee 2 


Column 


UL 219 


20 
31-40 
41-50 


SAI, 


DAoC 


k= 218, 


Format 


alee 
PiOe4 
re 


POs 
Pike 


> > 


ih 


F10.4 


Eon 4 


Tl 


F10.4 


ing 


Card 


Card 


8 


Description or Value 
5 
density lbs cu. fi. 


higher coefficients used as 
in Card 5 


2 
6 


time, seconds, at which the 
internal fluid temperature 
is that entered in columns 
21-30 


temperature of internal 
fluray oF... Type 6.-cards 
govern the time-temperature 
relationship of the internal 
fluid and as many as thirty 
may be used. These can be 
denoted as Cards 9b, 9c, 
etc. 


Cara He 


7 


Initial uniform temperature 
of pipe, °F 


Gard i 


8 Number eight (8) signals 
computer to begin computation 


Part II of the data deck must now be repeated using the 


same formats but this time using those parameters for the 


50 


connecting pipe (i.e., the pipe identified in Card 2). 
These cards follow directly after the card with the number 


eight (8) punched in column 1, and they must also end with 


paecard with 8 pinched in column I 


onl 


A(T) 

A3 
AA,BB3,RR 
AAA 

AENDX 
AENDMO 


AL (I) 


ARG 
ALPHA 
AYE 

B 

B3 
BBB3 
BBX 


BEST 


BETA 


BIG 


BUZZ (TI) 


APPENDIX B 


Glossary of Symbols Used in Program 


Listing of program starts on page 59 


cross sectional area at center number I 
inner radius of pipe 

quantities useful in calculating /Trdr 
quantity useful in calculating SIGMAR 
absolute value of ENDX 

absolute value of ENDMO 


Key (seo eh) =e NSiO(1)))) 7a 
useful combination 


value for which integral is evaluated 
coefficient of linear thermal expansion 
Inpue LEremecdara searea 

Input» £ Gcomeedeaneai-c 

outer radius of pipe 

quantity useful in calculating SIGMAR 
value of BETA times z 


maximum significant stress for a given 
time 





largest value of RT, TZ and ZR 


Variable used to store radius of max ?num 
Significant stress for each time interval 


a2 


BX 


C3 


cEeC 
D 


D3 


DENOM 
DENS (Z) 


DENA, DENB, 
DENC, DEND 


DEL 
DELTAL 
DELTA2 
BE TT. 
Dy Sd 


DONE, DTWO 


DD 


DDD 
DR 


ipa: 


ELAS 
EMINF 
ENDX 


ENDP 


BBX 
iMmput from dara Card 


Valbtew@ort ) Lrar from inner radius to some 
radius, R 


quantity useful in calculating SIGMAT 
input from data card 


Valwtewor frrdr fromminner radius to 
outer radius 


a quantity useful in calculating TAVG 
metal density at temperature z 


coefficients for determining DENS (Z) 


plate flexural rigidity = Flaeey sacar) 
AT, for USAS B31.7 Nuclear Power Piping Code 


AT. for USAS B31.7 Nuclear Power Piping Code 
absolute value of DELTAL 
absolute value of DELTA2 


quantities used in interpolating in 
subroutine FLUTEM 


quantity useful in determining radial 
displacement 


quantity useful in obtaining SIGMAZ 
thickness of a layer 

time interval 

input from data card 

modulus of elasticity 

bending moment remote from pipe junction 
bending moment as a function of Zz 


dislocation bending moment 


ayo 


ENDMO 

EPS2 

EF 

3 

FELTAL 

FPELTA2 

FLU (I )yaenwert ) 


FPLUDEERC Ae, ) 


FRAC 


Pou 
G 
GAKA 


GAUSS2 


GLUG1L 


GLUG2 


H3 


HO 


IDIMT 


IDIMX 


total bending moment at pipe junction 
unit elongation in axial direction 
input from data card 

quantity useful in calculating SIGMAR 
cumulatively largest value of DETT1 
cumulatively largest value of DETT2 
intertorymexterrouet | uid temperature 
subroutine entered with time Z and 
giving A,B as corresponding inside, 
outside fluid temperatures 

numerical constant used in making 
initial estimate of inner surface 
Cemperature 

interior fluid temperature 

quantity useful in calculating SIGMAT 
radius of maximum Significant stress 
Subroutine used for finding coefficients 
of polynomial for temperature distri- 
pwerton 

value of DETT1 at time TTl 

value of DETT2 at time TT2 

interior surface film coefficient 
quantity useful in calculating SIGMAZ 
external surface film coefficient (taken 
as zero corresponding to perfect 


LnSsiWeat | On 


number of coefficients in temperature 
polynomial 


number of coefficients in an associated 
temperature polynomial 


3 Reference to exterior fluid temperature is for a planed 
EUEUMe (Capability of —BEDIY!; NO such quantity 1s actually 
used in the calculations herein described. 
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UM 
JAP 
JUAY 


K (Z) 


KA,KB,KC,KD 
rea (T ) 
BORDER 

LL 

MO 


MOM 


MOP 


Nj Nj .N3,Ny 


NO 
NCOUNT 


NFLAG 


NPRINT 


NSTAR 
OD 


E 


ES 


EOLrS 
eon Ty 


PVAL 


IDIMX + l 
counter used in storing X3(I) 
integer parameter for preparing headings 


metal thermal conductivity at tempera- 
ture 2 


coefficients used in determining K(Z) 
K(X(1I))*A(I), useful combination 

degree of polynomial desired 

integer parameter for preparing headings 
LORDE Ret 1 


counter used in storing significant 
SEreSS fon Gacheerme interval 


counter used in storing significant 
stress of adjoining pipe for each time 
interval 


input from data card 


number of points used to make up poly- 
nomial for temperature distribution 


1 (in the present application of BETTY1) 
used to distinguish two joining pipes 


number of time intervals between lines 
of printed output 


option parameter in data deck 
outside diameter of pipe 


a useful parameter which takes care of 
conversion units 


quantity useful in obtaining end moment 
and radial force 


Poisson's ratio 
subroutine used to integrate polynomial 


subroutine used to evaluate /Trdr 


ao 


RAD (I) 


ROCK 


ne 


RY 


S3 


SHEAR 


SIG (I) 


SIGMAR 


5SIGMAT 


SIGMAZ 


Sire rx 


SIGZX 


SEGk 
SIGT 
SIGZ 


SPH AZ)) 


pene, SPH, 


BeHCy 2 oF HD 
SSIG 


Sit 


quantity useful in calculating end 
moment and radial shear force 


radial coordinate of pipe 


variable used to store radius of sig- 
bios Sigs 


radius of maximum significant stress for 
a given time 


absolute value of STl - ST2 


radial distance measured from mid-surface 
of pipe 


quantity useful in calculating end 
moment and radial force 


radial force aswa LuUNnctiOon Of axial 
coordinate 


variable used to store maximum signifi- 
cant stress for each time interval 


radial stress remote from pipe junction 


circumferential stress remote from pipe 
JUNC rem 


axial stress remote from junction 


circumferential stress due to bending 
moment and radial force 


axial stress due to bending moment and 
radial force 


total radial stress 

total circumferential stress 

total axial stress 

metal specific heat at temperature Z 
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maximum Significant stress encountered 


principal stress 1 
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SURFCO (Z) 
SURA, SURB, 
SuURC, SYRD 
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TIMA 
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TMAX 
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TVAL 


UAY (I) 
UU 


UUU 


VO 


W 


principal stress 2 
principal stress 3 


variable used to store significant 
Stress for a given radius 


inner surface film coefficient at 
temperature 2 


coefficients for calculating SURFCO(Z) 


problem-time 
shearing stress 


time at which maximum significant stress 
in encountered 


variable used to store time of maximum 
Significant stress for each time interval 


problem-end-time 
value of DT for first pipe 


time for which value of FELTAI1 is 
recorded 


time for which value "of 9FELTA2 is 
recorded 


temperature as obtained from polynomial 
for temperature distribution 


subroutine used to evaluate TT3 


radial displacement remote from pipe 
JUNGELOM 


alphabetic data for headings 

value of right hand side of equation (39) 
value of right hand side of equation (40) 
deflection at pipe junction 

radial force at pipe junction 


quantity useful in calculating end 
moment and radial force 
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wall thickness 


temperature at center I 


average temperature throughout thick- 
ness of pipe 


value of z where bending moment is 
extremal 


data point denoting radius used in 
fitting a polynomial 


coefficient for polynomial for tem- 
perature distribution 


outside fluid temperature 
newly calculated value of X(I) 


coefficients of polynomial for tem- 
perature after integration 


a useful quantity for calculating XAVG 
used in FLUTEM 


absolute value of ST3 - STl 


Recall that there are also a number of other variables with 


names similar to those above but with a terminal l. 
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